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ORIGINAL ARTICLE 

Oxidative stress induces mitochondrial dysfunction in a subset 
of autistic lymphoblastoid cell lines 

S Rose, RE Frye, J Slattery, R Wynne, M Tippett, S Melnyk and SJ James 

There is an increasing recognition that mitochondrial dysfunction is associated with autism spectrum disorders. However, little 
attention has been given to the etiology of mitochondrial dysfunction and how mitochondrial abnormalities might interact with 
other physiological disturbances such as oxidative stress. Reserve capacity is a measure of the ability of the mitochondria to 
respond to physiological stress. In this study, we demonstrate, for the first time, that lymphoblastoid cell lines (LCLs) derived from 
children with autistic disorder (AD) have an abnormal mitochondrial reserve capacity before and after exposure to reactive oxygen 
species (ROS). Ten (44%) of 22 AD LCLs exhibited abnormally high reserve capacity at baseline and a sharp depletion of reserve 
capacity when challenged with ROS. This depletion of reserve capacity was found to be directly related to an atypical simultaneous 
increase in both proton-leak respiration and adenosine triphosphate-linked respiration in response to increased ROS in this AD LCL 
subgroup. In this AD LCL subgroup, 48-hour pretreatment with ^-acetylcysteine, a glutathione precursor, prevented these 
abnormalities and improved glutathione metabolism, suggesting a role for altered glutathione metabolism associated with this 
type of mitochondrial dysfunction. The results of this study suggest that a significant subgroup of AD children may have alterations 
in mitochondrial function, which could render them more vulnerable to a pro-oxidant microenvironment as well as intrinsic and 
extrinsic sources of ROS such as immune activation and pro-oxidant environmental toxins. These findings are consistent with the 
notion that AD is caused by a combination of genetic and environmental factors. 
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INTRODUCTION 

The autism spectrum disorders (ASDs) are a heterogeneous group 
of neurodevelopmental disorders defined by impairments in 
communication and social interaction along with restrictive and 
repetitive behaviors.^ An estimated 1 out of 88 individuals in the 
United States are currently affected with an ASD and the incidence 
continues to rise.^ 

Mitochondrial dysfunction has become increasingly recognized 
as a major physiological disturbance in ASD.^ However, the 
etiology of mitochondrial dysfunction is not known. Indeed, 
although mitochondrial deoxyribonucleic acid mutations are 
commonly found in classical mitochondrial disease (MD), such 
mutations are only found in 23% of ASD children diagnosed with 
MD.^ This raises the possibility of acquired mitochondrial 
dysfunction as mitochondrial damage can result from environ- 
mental exposures implicated in ASD such as heavy metals,"^"^ 
exhaust fumes,^ polychlorinated biphenyls^ or pesticides.^°'" 
Alternatively, mitochondria can be damaged by endogenous 
stressors associated with ASD such as elevated proinflammatory 
cytokines resulting from an activated immune system^^"^"^ or 
other conditions associated with oxidative stress.^ ^'^^ The notion 
of an acquired mitochondrial disorder is supported by a recent 
twin study, which concluded that the environment contributes a 
greater percent of the risk of developing autistic disorder (AD, 
55%) as compared with genetic factors (37%) with these factors 
contributing about equally for the broader ASD diagnosis.^^ 

Oxidative stress may be a key link between mitochondrial 
dysfunction and ASD as reactive oxygen species (ROS) generated 



from pro-oxidant environmental toxicants^"" and activated 
immune cells^'^^ can result in mitochondrial dysfunction.^ Four 
independent case-control studies have documented oxidative 
stress and oxidative damage in plasma, immune cells and post- 
mortem brain from ASD children. ^^"^^ Interestingly, resting 
peripheral blood mononuclear cells and activated monocytes 
derived from children with ASD demonstrate a significant 
decrease in glutathione redox balance reflecting an intracellular 
deficit in glutathione-mediated antioxidant and detoxification 
capacity.^^ 

Lymphoblastoid cell lines (LCLs) are readily available in national 
and international biorepositories for many psychiatry and 
neurological diseases, including autism, and have been used as 
biological models for psychiatric and neurological disease, 
including the examination of genetic and metabolic aspects of 
these diseases. Developing a model of mitochondrial function for 
children with ASD using these cell lines would allow the detailed 
study of mitochondrial function as well as its associated metabolic 
abnormalities. We have previously demonstrated the cellular and 
mitochondrial redox imbalances in LCLs derived from children 
with AD.^° Thus, in this study we hypothesized that LCLs derived 
from patients with AD are vulnerable to ROS, such that excessive 
intracellular ROS results in mitochondrial dysfunction. Further- 
more, we hypothesize that this dysfunction only affects a subset of 
LCLs derived from children with AD. To this end, we examined 
mitochondrial respiratory activity in LCLs derived from AD children 
and unaffected control individuals. We demonstrate atypical 
changes in mitochondrial respiration in AD LCLs when exposed 



Department of Pediatrics, Arkansas Children's Hospital Research Institute, University of Arkansas for Medical Sciences, Little Rock, AR, USA. Correspondence: Dr RE Frye, 
Department of Pediatrics, Arkansas Children's Hospital Research Institute, University of Arkansas for Medical Sciences, Slot 51 2-41 B, Room R4041, 13 Children's Way, Little Rock, 
AR 72202, USA. 
E-mail: REFrye@uams.edu 

Received 6 November 2013; revised 31 January 2014; accepted 2 February 2014 



Mitochondrial dysfunction in autism cell lines 

S Rose et al 



to ROS in a subgroup of AD LCLs. We reveal the significance of 
glutathione metabolism in this effect by showing that pretreat- 
ment with ^-acetylcysteine (NAC), a glutathione precursor, 
prevents mitochondrial vulnerability in the AD LCL subgroup that 
exhibit atypical responses to increased ROS. 



MATERIALS AND METHODS 

Lymphoblastoid cell lines and culture conditions 
Twenty-two LCLs derived from white males diagnosed with AD chosen 
from pedigrees with at least one affected male sibling (mean/s.d. age 
7.8 ±3.1 year) were obtained from the Autism Genetic Resource Exchange 
(Los Angeles, CA, USA) or the National Institutes of Mental Health 
(Bethesda, MD, USA) center for collaborative genomic studies on mental 
disorders. Fourteen control LCLs derived from healthy white male donors 
with no documented behavioral or neurological disorder (mean/s.d. age 
27.7 ±9.1) were obtained from Coriell Cell Repository (Camden, NJ, USA) or 
the National Institutes of Mental Health. AD LCLs and control LCLs were 
randomly paired with each other (Supplementary Table SI). Young adult 
controls are standard in LCL studies due to low availability of younger 
ages.^° On average, cells were studied at passage 12, with a maximum of 
15. Genomic stability is very high at this low passage.^^'^^ Cells were 
maintained in RPMI 1640 culture medium with 15% FBS and 1% penicillin/ 
streptomycin in a humidified incubator at 37 °C with 5% CO2. 

Seahorse assay 

To measure mitochondrial function, we used the state-of-the-art Seahorse 
Extracellular Flux (XF) 96 Analyzer (Seahorse Bioscience, North Billerica, MA, 
USA). On the day of the assay, XF-PS plates were coated with 13 |al poly-D- 
lysine (50 |ag mT^) for 2 h and washed twice with cell culture-grade water. 
One hour before the assay, cells were seeded onto coated 96-well XF-PS 
plates at a density of 1.1 x 10^ cells/well in Dulbecco's modified Eagle's 
medium XF assay media (unbuffered Dulbecco's modified Eagle's medium 
supplemented with 1 1 mivi glucose, 2 mivi L-glutamax and 1 mivi sodium 
pyruvate). Replicate/well varied from 4 to 8 depending on the number of 
cells available. Titrations determine the optimal concentrations of 
oligomycin (1.0|aivi), carbonyl cyanide-p-trifluoromethoxyphenyl-hydrazon 
(0.3 hm), antimycin A (0.3 |aivi) and rotenone (1.0 |aivi). By sequentially adding 
these pharmacological agents to the respiring cells, we measured the basal 
respiration, adenosine-5'-triphosphate-linked respiration, PLR and RC (see 
Figure 1). 

Redox challenge 

ROS was increased in vitro by exposing cells to increasing concentrations 
of the redox cycling agent DMNQ for 1 h before the Seahorse assay. DMNQ 
enters cells and generates both superoxide and hydrogen peroxide similar 
to that generated by NADPH oxidase in vivo}^ A 5mgmr^ DMNQ 
solution was diluted in Dulbecco's modified Eagle's medium XF assay 
media into 10X stocks and added to cells in an XF-PS plate and incubated 
for 1 h at 37 °C in a non-C02 incubator. The concentrations of DMNQ were 
optimized as 5, 10, 12.5 and 15 [jm. 

NAC rescue 

To determine whether pretreatment with a glutathione precursor could 
rescue atypical responses to the ROS challenge, AD LCLs were plated in 
T25 flasks at a density of 5 x 1 0^ cells ml" ^ in culture media with or without 
1 mM NAC for 48 h. Cells were washed twice in Dulbecco's modified Eagle's 
medium XF media to remove any remaining NAC before DMNQ treatment 
and Seahorse assays. The control LCLs were not treated with NAC as they 
demonstrate normal glutathione metabolism. 

Glutathione measurements 

Approximately 2x10^ viable cells were pelleted and snap-frozen on dry 
ice. Samples were stored at - 80 °C until HPLC quantification of intracellular 
free GSH and GSSG.^^ Thawed cells were lysed by 3 s sonication in 1 12.5 |al 
ice-cold phosphate-buffered saline followed by the addition of 37.5 |al ice- 
cold 10% meta-phosphoric acid. This mixture was incubated for 30 min on 
ice followed by centrifuging for 15 min at 18 000xg at 4°C. Results are 
expressed as nmolmg"^ protein using BCA Protein Assay Kit (Pierce, 
Rockford, IL, USA). 
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Figure 1. The Seahorse assay. Oxygen consunnption rate is nneasured 
before and after adding pharnnacological agents to respiring cells. 
Measurennent of oxygen consunnption over 6 nnin is nnade repeat- 
edly. Three nneasurennents are nnade and averaged to provide 
reliable nneasurennents. For the first 18 nnin, total cellular oxygen 
consunnption is nneasured. Basal respiration can be calculated fronn 
this quantity by subtracting non-nnitochondrial respiration. Next 
oligonnycin, an inhibitor of adenosine-5'-triphosphate (ATP) respira- 
tion, is added and a nneasurennent of this is nnade over the next 
18 nnin. This quantity can be subtracted fronn the total cellular 
oxygen consunnption to deternnine ATP-linked respiration and non- 
nnitochondrial respiration can be subtracted fronn this quantity to 
obtain proton-leak respiration. Next carbonyl cyanide-p-trifluoro- 
nnethoxyphenyl-hydrazon, a protonophore, is added and a nneasure- 
nnent of this is nnade over the next 18 nnin. The protonophore 
collapses the inner nnennbrane gradient by nnaking the inner 
nnennbrane pernneable to protons. This drives the electron transport 
chain to function at its nnaxinnunn rate. Subtracting non- 
nnitochondrial respiration fronn this quantity produces a nneasure 
of nnaxinnunn respiratory capacity. Finally antinnycin A, a connplex III 
inhibitor, and rotenone, a connplex I inhibitor, are added to shut 
down electron transport chain function. The resulting nneasurennent 
over the next 18 nnin represents non-nnitochondrial respiration, a 
nneasurennent that can be used with the other nneasurennents to 
calculate respiratory paranneters. Finally, reserve capacity is calcu- 
lated by subtracting basal respiration fronn nnaxinnunn respiratory 
capacity. 



Analytic approach 

A mixed-effects regression^^ was conducted via SAS version 9.3 (Cary, NC, 
USA) 'gimmix' procedure. The mixed-effects models allowed data from 
each AD LCL to be compared with the paired control LCL run on the same 
plate. The mitochondrial respiratory measurement was the response 
variable with a between-group dichotomous effect (AD vs control, AD 
-i-NAC vs control) and within-group repeated factors of DMNQ concentra- 
tion (modeled as a multilevel factor) as well as the interaction between 
these effects. We present the overall difference between the two 
comparison groups (group effect), the overall effect of the DMNQ 
concentration (DMNQ effect), and whether the effect of DMNQ concentra- 
tion was different between the two groups (DMNQ x group interaction). 
Random effects included the intercept and DMNQ. F-tests were used to 
evaluate significance. Planned post hoc orthogonal contrasts were used 
when the interaction was significant. 

Differences in glutathione measurements between control and AD LCLs 
without DMNQ exposure were analyzed using a similar mixed-effect 
regression model. A general linear model was used to verify the DMNQ 
effect on AD and control LCLs as these LCLs were not matched. DMNQ was 
treated as a continuous variable since a dose response effect was 
expected. 

Cluster analysis was conducted using Ward's technique.^^ Ward's 
technique defines the distance between clusters in terms of the between 
cluster variability to the within cluster variability. By examining the 
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dendogram and several statistics (pseudo-F and t), a judgment is made 
about the number of clusters.^° 

RESULTS 

Verification of 2,3-climethoxy-1,4-napthoquinone (DMNQ) effect 
on LCL glutathione redox state 

Glutathione concentrations were measured in three control and 
five AD LCLs at 2,3-climethoxy-1,4-napthoquinone (DMNQ) con- 
centrations of 0, 1, 5, 10, 12.5 and 15 mm. DMNQ significantly 
decreased reduced glutathione (GSH) (F(1, 35) -52.45, P< 0.001) 
and GSH/oxidized glutathione (GSSG) (F(1, 35) = 30.21, P< 0.001) 
and increased GSSG (F(1, 35) = 13.80, P < 0.001) in a linear fashion 
(see Supplementary Figure SI). These changes were not different 
across groups. 

Mitochondrial function in AD LCLs with ROS challenge 
Basal respiration was greater overall in AD LCLs (F(1,1 05) = 10.02, 
P< 0.005) and significantly increased as DMNQ increased (F 
(4,84) = 6.22, P < 0.0005) (Figure 2a). This increase was significantly 
greater in AD LCLs (F(4,105) = 5.20, P< 0.001), due to a 
significantly higher basal respiration in AD LCLs at 10 hm (f 
(105) = 4.42, P< 0.0001) and 12.5 mm (f(105) = 3.21, P< 0.005). 
Basal respiration is composed of proton-leak respiration (PLR) 
and Adenosine-5'-triphosphate-Linked Respiration (ALR), which 
we examined separately. 

Overall ALR was higher for AD LCLs (F(1,105) = 5.84, P<0.05) 
(Figure 2b). ALR increased as DMNQ was increased ( F (4,84) = 5.44, 
P< 0.001) with this increase significantly greater in AD LCLs (F 
(4,1 05) = 3.97, P< 0.005), due to significantly higher ALR in AD 
LCLs at 10 mm (f(105) = 3.86, P< 0.0005) with borderline signifi- 
cance at 1 2.5 MM (f(1 05) = 1 .78, P < 0.1 0). 

PLR was overall greater in AD LCLs (F(1, 105) = 21.05, P< 0.0001) 
(Figure 2c) and significantly increased as DMNQ increased (F 
(4,84) = 11.31, P< 0.0001). This increase was significantly greater 
for AD LCLs (F(4,1 05) = 4.29, P = 0.005), due to a significantly higher 
PLR in AD LCLs at 10 mm (f(105) = 3.75, P< 0.0005) and 12.5 mm (f 
(1 05) = 4.71, P< 0.0001) (Figure 2c). 



Maximal respiratory capacity (MRC) was overall significantly 
higher in AD LCLs (F(1,105) = 6.76, P = 0.01). Both LCL groups 
showed a significant (and similar) decrease in MRC as DMNQ 
increased ( F (4,84) = 35.88, P< 0.0001) (Figure 2d). 

As DMNQ increased, reserve capacity (RC) decreased (F 
(4,84) = 50.52, P< 0.0001) with this decrease significantly greater 
for AD LCLs (F(4,105) = 6.16, P< 0.0005). RC of AD LCLs was 
significantly higher than control LCLs without DMNQ (f 
(105) = 2.86, P = 0.005) and at 5 mm (f(105) = 2.73, P<0.01) but 
then sharply dropped to be significantly lower than control LCLs 
at 10 mm DMNQ (f(105) = 2.44, P<0.05) (Figure 2e). 



Defining subgroups of AD LCLs 

As AD and control LCLs differed markedly in the changes in ALR 
and PLR with DMNQ challenge, we examined whether the 
changes in these respiratory parameters could differentiate AD 
LCL subgroups. As the increase in ALR and PLR peaked at 10 mm 
DMNQ, the slope of the change in ALR and PLR from 0 to 10 mm 
DMNQ was calculated and entered into a cluster analysis. The 
cluster analysis divided the LCLs into two groups: AD-N (n = 12) 
and AD-A (n = 10) (pseudo-F 17.1, pseudo t-squared 18.3) (see 
Figure 3). The dendogram (not shown) demonstrated clear 
differences between these groups. 

The most striking difference between the two groups was the 
relationship between the change in PLR and ALR. For the AD-N 
group, there was a significant negative relationship between PLR 
and ALR (Figure 3a, r=-0.77, P<0.01, green line) such that an 
increase in PLR was associated with a decrease in ALR (and vice 
versa). If the two AD-N outliers were removed, the correlation 
remained high for AD-N LCLs (r = -0.86, P<0.01; blue line). 
Overall, there was little change in PLR and ALR with the DMNQ 
challenge in the AD-N subgroup (Figures 3b and c, green dashed 
lines). 

For the AD-A subgroup, the relationship between PLR and ALR 
was different than the AD-N subgroup. Specifically, PLR and ALR 
both increased together (Figure 4a, r = 0.44, P=NS, red line; 
Figures 3d and e, red dashed lines), rather than having an inverse 
relationship. 
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Figure 2. Lymphoblastoid cell lines (LCLs) derived from children with autistic disorder (AD) demonstrate differences in mitochondrial function 
as compared with control LCLs at baseline and after exposure to the redox cycling agent 2,3-dimethoxy-1,4-napthoquinone (DMNQ) at four 
concentrations (5, 10, 12.5 and 15 [im) one hour before the assay, (a) Basal respiration increases as DMNQ concentration increases in the AD 
LCLs and becomes significantly higher in the AD LCLs at 10 and 12.5 \iM DMNQ; (b) adenosine-5'-triphosphate-linked respiration increases as 
DMNQ concentration increases in the AD LCLs and becomes significantly higher in the AD LCLs at 10|iM DMNQ; (c) proton-leak respiration 
increases as DMNQ concentration increases in the AD LCLs and becomes significantly higher in the AD LCLs at 10 and 12.5|iM DMNQ; (d) 
maximum respiratory capacity decreased as DMNQ increased for both AD and control LCLs but overall AD LCLs demonstrated a higher 
maximum respiratory capacity; (e) reserve capacity decreases as DMNQ increases for both AD and control LCLs but the decline in reserve 
capacity is much sharper for the AD LCLs as compared with the control LCLs due to the fact that reserve capacity is significantly higher in the 
AD LCLs at low DMNQ concentrations but becomes significantly lower in the AD LCL at higher DMNQ concentrations. 
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Figure 3. Clustering of lymphoblastoid cell lines (LCLs) derived from children with autistic disorder (AD) into two subgroups: the AD-N 
subgroup has mitochondrial respiratory profiles similar to controls and the AD-A subgroup has atypical mitochondrial respiratory profiles, (a) 
AD-N cases (green circles) and AD-A cases (red diamonds) represent the two subgroups. The AD-N subgroup demonstrates a negative 
correlation between adenosine-5'-triphosphate (ATP)-linked respiration and proton-leak respiration (r = -0.77, P<0.01, green line) whereas 
the AD-A group demonstrates a positive correlation between ATP-linked respiration and proton-leak respiration (r = 0.44, P=NS, red line). If 
the two outliers are removed from the AD-N group, the correlation is still significant (r=-0.86, P < 0.01; blue line), (b) Individual (thin lines) 
and overall (thick green dashed line) change in ATP-linked respiration for the AD-N groups. Notice that there is little change overall, (c) 
Individual (thin lines) and overall (thick green dashed line) change in proton-leak respiration for the AD-N groups. Notice that there is little 
change overall, (d) Individual (thin lines) and overall (thick red dashed line) change in ATP-linked respiration for the AD-A groups. Notice that, 
overall, ATP-linked respiration increases with DMNQ concentration, (e) Individual (thin lines) and overall (thick red dashed line) change in 
proton-leak respiration for the AD-A groups. Notice that, overall, proton-leak respiration increases with DMNQ concentration. DMNQ, 2,3- 
dimethoxy-1,4-napthoquinone; NS, not significant. 



Q 
< 



ATP-Linked Respiration 

a 

^80 



Q. 



-♦-Control 

-AD-N 



0 5 10 15 
DMNQ Concentration (uM) 




< 
I 

Q 
< 



-♦-Control * p<0.05 
-^AD-A T p<0.01 



0 5 10 15 
DMNQ Concentration (uM) 



Proton Leak Respiration Maximum Respiratory Capacity 



J 40 



Reserve Capacity 



1 1"^-+^ t 



20 



o 



^ 40 

■5 
E 

f 20 

9) 





-♦-Control 

-•-AD-N 








1 1 1 1 
0 5 10 15 
DMNQ Concentration (uM) 


V 

-♦-Control * T 


1 r 


T p<0.01 
t p<0.0001 




0 5 10 15 
DMNQ Concentration (uM) 



0 5 10 15 
DMNQ Concentration (uM) 



0 5 10 15 
DMNQ Concentration (uM) 




0 5 10 15 
DMNQ Concentration (uM) 



0 5 10 15 
DMNQ Concentration (uM) 



Figure 4. Seahorse respiratory measurements in two subgroups of lymphoblastoid cell lines (LCLs) derived from children with autistic disorder 
(AD) as compared with control LCLs at baseline and after exposure to the redox cycling agent 2,3-dimethoxy-1,4-napthoquinone (DMNQ) at 
four concentrations (5, 10, 12.5 and 15 |iM) one hour before the assay. Overall, the AD-A subgroup (e-h) parallels the differences between the 
AD and control LCLs found in the overall analysis whereas the AD-N subgroup (a-d) demonstrates similar mitochondrial response between 
the AD LCLs and control LCLs. For the AD-N subgroup (a) adenosine-5'-triphosphate(ATP)-linked respiration, (b) proton-leak respiration, (c) 
maximum respiratory capacity and (d) reserve capacity are similar between the AD and control LCLs. For the AD-A subgroup, (e) ATP-linked 
respiration increases as DMNQ concentration increases in the AD-A LCLs and becomes significantly higher in the AD-A LCLs at IGjiivi DMNQ; 
(f) proton-leak respiration increases as DMNQ concentration increases in the AD-A LCLs and becomes significantly higher in the AD-A LCLs at 
10 |iM and 12.5 |iM DMNQ; (g) maximum respiratory capacity decreased as DMNQ increased for both AD-A and control LCLs but overall AD-A 
LCLs demonstrated a higher maximum respiratory capacity; (h) the decline in reserve capacity is much sharper for the AD-A LCLs as compared 
with the control LCLs due to the fact that reserve capacity is significantly higher in the AD-A LCLs at low DMNQ concentrations but becomes 
significantly lower in the AD-A LCL at higher DMNQ concentrations. 



As a comparison, the relationship between the PLR and ALR was 
examined for the control LCLs. The relationship between PLR and 
ALR for the control LCLs was found to be significantly positive 
(r = 0.79, P<0.01; not shown). 



Mitochondrial function in AD LCLs subgroups with ROS challenge 
To better understand the differences between the AD LCL 
subgroups, we compared the AD LCLs with their paired control 
LCLs within each subgroup. Basal respiration is not shown as it is 
the combination of ALR and PLR. 
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Figure 5. Seahorse respiratory nneasurennents in two subgroups of lynnphoblastoid cell lines (LCLs) derived fronn children with autistic disorder 
(AD) after 48 h incubation with 1 nmvi of A/-acetyl-cysteine (NAC) as connpared with control LCLs at baseline and after exposure to the redox 
cycling agent 2,3-dinnethoxy-1,4-napthoquinone (DMNQ) at four concentrations (5, 10, 12.5 and 15|iM) one hour before the assay. Overall, 
atypical changes in nnitochondrial function with increased DMNQ seen in the AD-A subgroup (e-h) are rescued with NAC treatnnent whereas 
NAC does not alter the dynannics of nnitochondrial function in the AD-N subgroup (a-d). For the AD-N subgroup (a) adenosine-5'-triphosphate 
(atp)-linked respiration, (b) proton-leak respiration, (c) nnaxinnunn respiratory capacity and (d) reserve capacity are sinnilar between the AD and 
control LCLs. For the AD-A subgroup, changes (e) ATP-linked respiration, (f) proton-leak respiration, (g) nnaxinnunn respiratory capacity and (h) 
reserve capacity with increasing DMNQ did not differ between the AD-A and control LCLs. However, (e) ATP-linked respiration, (f) proton-leak 
respiration and (g) nnaxinnunn respiratory capacity were found to be overall greater in the AD-A LCLs as connpared with control LCLs. This 
change in the overall function of the AD-A LCLs with NAC pretreatnnent nornnalized the reserve capacity differences between the AD-A and 
control LCLs. 



AD-N LCLs. ALR, PLR, MRC and RC significantly changed as DMNQ 
increased in both AD-N and control LCLs (F(4,40) = 3.22, P < 0.05; F 
(4,40) = 3.1 3, P<0.05; F(4,40) = 15.92, P< 0.0001; F(4,40) = 22.74, 
P< 0.0001, respectively) (Figures 4a-d), but this change was not 
significantly different between the LCL groups. The AD-N LCLs did 
not dennonstrate any overall differences from controls for any of 
these respiration parameters. 

AD-A LCLs. Overall, ALR was higher for AD-A LCLs (F(1 ,50) = 6.46, 
P = 0.01) (Figure 4e). ALR significantly increased as DMNQ 
increased in both LCL groups (F(4,40) = 3.37, P<0.05) with this 
increase significantly greater for AD-A LCLs (F(4,50) = 2.92, 
P<0.05), due to a greater ALR in the AD-A LCLs at 10 mm (f 
(50) = 3.23, P < 0.01 ) and 1 2.5 mm (f(50) = 2.02, P < 0.05). 

PLR was overall higher for the AD-A LCLs (F(1, 50) = 24.31, 
P< 0.0001) (Figure 4f). Increasing DMNQ resulted in a significant 
increase in PLR ( F (4,40) = 1 2.52, P< 0.0001) with this increase 
significantly greater for the AD-A LCLs (F(4,50) = 3.60, P = 0.01 ), due 
to PLR being significantly higher in the AD-A LCLs at 10 mm (f 
(50) = 4.48, P < 0.0001 ) and 1 2.5 mm (f(50) = 3.56, P < 0.001 ). 

MR decreased as DMNQ increased (F(4,40) = 1 9.63, P< 0.0001) 
and was overall significantly greater for the AD-A LCLs (F 
(1,50) = 9.58, P< 0.005) (Figure 4g). 

RC decreased as DMNQ increased (F(4,40) = 28.41, P< 0.0001) 
with this decrease significantly more marked for AD-A LCLs (F 
(4,50) = 6.08, P = 0.0005) (Figure 4h). This was due to the RC being 
significantly higher in AD-A LCLs at low DMNQ concentrations 
(Omm f(50) = 3.18, P< 0.005; 5 mm f(50) = 2.56, P = 0.01) and then 
plummeting so that RC was significantly lower in the AD-A LCLs at 
1 0 UM (1 0 MM f(50) = 2.1 6, P < 0.05). 



A/-acetyl-l-cysteine rescues AD LCL RC abnormalities 

AD LCLs were pre-treated with NAC to determine whether 

improving glutathione metabolism could normalize mitochondrial 

function. 

AD-N NAC-treated LCLs. ALR, PLR, MRC and RC significantly 
changed as DMNQ increased ( F (4,40) = 3.04, P<0.05; F 
(4,40) = 3.1 3, P<0.05; F (4,40) = 15.92, P< 0.0001; F (4,40) = 22.74, 
P< 0.0001, respectively) (Figures 5a-d), but this change was not 
significantly different between the LCL groups. The AD-N LCLs did 
not demonstrate any overall differences from their matched 
controls for any respiration parameter. 

AD-A NAC-treated LCLs. ALR, PLR, MRC and RC significantly 
changed as DMNQ increased ( F (4,40) = 4.09, P<0.01; F 
(4,40) = 3.1 9, P<0.05; F(4,40) = 15.28, P< 0.0001; F (4,40) = 21 .67, 
P< 0.0001, respectively) (Figures 5e-h), but, this change was not 
significantly different between the AD-A and control LCLs. Instead, 
ALR, PLR, MRC were overall significantly higher in AD-A LCLs (F 
(1 ,50) = 1 0.1 9, P < 0.0001 ; F(1 ,50) = 1 5.08, P < 0.0005; F(1 ,50) = 4.25, 
P < 0.05, respectively). 



Intracellular glutathione metabolism in LCLs 
AD LCLs demonstrated lower GSH (F(1, 89) = 121.50, P< 0.001) and 
GSH/GSSG ratio (F(1, 89) = 234.60, P< 0.0001) and higher GSSG (F 
(1,89) = 420.11, P< 0.0001) as compared with control LCLs (see 
Supplementary Figure S2). 

Compared with control LCLs, the two AD LCL subgroups 
demonstrated significantly lower GSH (AD-A: F(1, 53) = 200.18, 
P< 0.0001; AD-N: F(1, 35) = 14.46, P< 0.001) and GSH/GSSG (AD- 
A: F(1,53) = 443.66, P< 0.0001; AD-N: F(1, 35) = 92.89, P< 0.0001) 
and significantly higher GSSG (AD-A: F(1, 53) = 1 14.19, P< 0.0001; 
AD-N: F(1,35) = 131.78, P< 0.0001), however neither GSH, GSSG 
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nor GSH/GSSG ratio were significantly different between the two 
AD LCL subgroups (see Supplementary Figure S2). 

Pretreatment with NAC increased intracellular GSH and the 
GSH/GSSG ratio and reduced GSSG in both the AD-A (GSH: F 
(1,53) = 783.69, P< 0.0001; GSSG: F(1, 53) = 810.73, P< 0.0001; 
GSH/GSSG: F(1, 53) = 1264.48, P< 0.0001) and AD-N (GSH: F 
(1,35) = 87.34, P< 0.0001; GSSG: F(1, 35) = 22.50, P< 0.0001; GSH/ 
GSSG: F(1,35) = 110.22, P< 0.0001) LCLs as compared with these 
LCLs before pretreatment (see Supplementary Figure S2). 



DISCUSSION 

This study examined the mitochondrial function in immune cells 
derived from children with AD at baseline and after exposure to 
an agent that increased ROS levels in vitro. Here, for the first time, 
we show that LCLs derived from children with AD exhibit 
significant abnormalities in mitochondrial respiration after expo- 
sure to increasing levels of ROS. Specifically, we demonstrate 
atypical simultaneous increases in ALR and PLR along with a sharp 
drop in RC when exposed to ROS in the AD LCLs as compared with 
the control LCLs. We then further demonstrate that these atypical 
responses were driven by a subset that comprised 44% of the AD 
LCLs. This subgroup also demonstrated an increase in MRC, further 
resulting in an increase in RC when the LCLs were not exposed to 
ROS. The atypical response to increased ROS in this AD LCL 
subgroup was prevented by an NAC pretreatment. Overall, this 
study demonstrates that at least a subset of LCLs derived from 
children with AD may have physiological abnormalities in 
mitochondrial function that are associated with abnormalities in 
oxidative stress. This evidence provides important insight into the 
potential pathophysiological mechanisms associated with AD and 
potential pathways for treatment. 

Differences in RC depletion in AD LCL subgroups with ROS 
exposure 

The adaptive response to ROS were divided into normal (AD-N) 
and abnormal (AD-A) subgroups. Exposure to ROS resulted in a 
more precipitous decrease in RC in the AD-A as compared with the 
AD-N LCLs. This is significant since reduced RC is linked to several 
diseases such as aging,^^ heart disease^^ and neurodegenerative 
disorders.^^'^^ In addition, when RC is depleted, apoptosis ensues 
resulting in reduced cell viability.^^ Differences in the dynamics 
between ALR and PLR in the AD LCL subgroups accounted for 
subgroup-specific responses. 

In the AD-N subgroup, neither ALR nor PLR changed 
significantly, on average, as ROS increased. On an individual LCL 
level, we found a negative relationship between the change in 
ALR and PLR with increasing ROS. AD-N LCLs demonstrated an 
increase in ALR or PLR in response to increasing ROS, whereas the 
respiratory parameter that did not increase (ALR or PLR) 
decreased. This resulted in stable basal respiration as basal 
respiration is the combination of ALR and PLR. This, in turn, limited 
the decrease in RC as RC is the difference between MRC and basal 
respiration. These observations suggest that AD-N LCLs respond to 
increased ROS with one of two mechanisms: either increasing PLR 
or ALR, but not both. 

For AD-A LCLs, the simultaneous increase in both ALR and PLR, 
on average and possibly on an individual level, resulted in a basal 
respiration increase which, when subtracted from MRC, results in a 
decrease in RC. Thus, this accounts for the greater decrease in RC 
in AD-A LCLs as compared with AD-N LCLs. This also suggests that 
ALR and PLR are required to work together to respond to 
increased ROS in the AD-A LCLs. Interestingly, the control LCLs 
demonstrated a significant positive relationship between ALR and 
PLR, yet did not demonstrate the same depletion in RC as the AD- 
A LCLs. This suggests that the AD-A are using the same 
mechanisms as the control LCLs to respond to the increases 
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ROS but cannot sustain this response properly, most likely 
because these ROS control mechanism are inadequate in the 
context of chronic elevations in ROS. Alternatively, it appears that 
the AD-N LCLs have adapted to control for acute increases in ROS 
differently than control, probably as a consequence of chronic 
elevations in ROS. 

RC depletion appears to be maximal at about 10|jm DMNQ. In 
fact ALR and PLR appear to peak at 10 mm DMNQ and then 
decreases at higher DMNQ concentrations. This decrease in ALR 
and PLR at high DMNQ concentrations most likely represents a 
failure of the mitochondria to function after RC is depleted. As 
depletion in RC results in reduced viability, it is likely that an 
apoptotic cascade is initiated when RC is depleted. This depletion 
most likely occurs faster as DMNQ concentration increases. As the 
LCLs are exposed to DMNQ for the same period of time regardless 
of the DMNQ concentration, it is likely that the LCLs exposed to 
DMNQ concentrations higher than 10 |jm will have greater cellular 
damage after the 1h exposure and this will be reflected in the 
measurements of the respiratory parameters during the assay. 

NAC rescues the atypical mitochondrial respiratory response in 
AD-A LCLs 

Pretreatment of the AD-A LCL subgroup with NAC improved 
glutathione redox status and improved mitochondrial respiration. 
Interestingly, NAC did not normalize respiratory parameters but 
rather increased the ability of the cell to produce adenosine 
triphosphate (ATP) and maintain an adequate maximum respira- 
tory capacity in the context of increased ROS. This suggests that 
an insufficient antioxidant capacity could contribute to the 
vulnerability of the mitochondrial to ROS and suggests that the 
AD-A LCL subset we have identified are particularly dependent on 
the cellular redox status. Interestingly, NAC may be efficacious in 
treating drug dependence,^^ schizophrenia,^^ bipolar 
depression,^^ trichotillomania^^ and may have a role in the 
treatment of Alzheimer's disease.^°'^^ Recent studies suggest that 
NAC has a role in supporting mitochondrial metabolism,^^ 
including mitigating oxidative stress-induced mitochondrial 
dysfunction,"^^'"^^ attenuating mitochondrial-related oxidative 
stress"^^ and restoring cognitive deficits in a mouse model with 
complex I deficiency.^^ In a recent small double-blind placebo- 
controlled trial, NAC significantly improved irritability and social 
cognition in ASD children."^^ 

Molecular mechanisms associated with the increase in ATP-linked 
respiration 

Elevations in ALR in the AD-A LCL subgroup is consistent with 
clinical reports of electron transport chain (ETC) overactivity in 
ASD children. Frye and Naviaux'^'' reported five ASD/MD children 
with complex IV overactivity and Graf et at^ reported a ASD/MD 
child with complex I overactivity. The fact that ALR increases in 
response to increased ROS, suggests that increased ATP produc- 
tion may be an important cytoprotective mechanism against ROS. 
MRC is a measure of the maximum ability of the ETC to generate 
ATP. Higher MRC in the AD-A subgroup is consistent with an 
overall increase in ALR and, again, suggests an overactivity of the 
ETC in the AD-A LCLs. 

Molecular mechanisms associated with the increased PLR 
Proton leak reduces the mitochondrial membrane potential which, 
in turn, decreases ETC ROS generation."^^ Proton leak is modulated 
by several mechanisms, including the adenine translocator and, in 
lymphocytes, uncoupling protein 2.^° Given that uncoupling 
protein 2 is upregulated by chronic oxidative stress^^'" and that 
AD LCLs have chronic elevations in ROS,^° it is possible that an 
increase in uncoupling protein 2 could be associated with the 
increase in PLR in the AD-A subgroup. However, as an increase in 
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mitochondrial membrane potential can increase both PLR and 
ALR,^^ an increase in mitochondrial membrane potential, poten- 
tially driven by ETC complex overactivity, could also account for 
the findings in the AD-A subgroup. Further research will be 
needed to clarify the molecular mechanisms leading to our 
findings. 

MD and dysfunction in ASD 

The nature and prevalence of MD in ASD is still under 
investigation. A recent meta-analysis found that 5% of children 
with ASD meet criteria for a classic MD but that 30+% of children 
in the general ASD population exhibit biomarkers consistent with 
MD.^ Recently, Frye^"^ demonstrated that 50+% of ASD children 
have biomarkers of MD that are consistently abnormal (that is, 
repeatable) and valid (that is, correlate with other MD biomarkers). 
In another study, 80% of the children with ASD demonstrated 
abnormal lymphocytes ETC function.^^ Interestingly some children 
with ASD/MD have ETC overactivity rather than ETC 
deficiencies"^^'"^^ and many ASD/MD cases do not manifest lactate 
elevation,'^'''^^'^^ a key biomarker commonly used to identify 
individuals with classic MD. This has raised the idea that children 
with ASD might have a type of mitochondrial dysfunction that is 
more prevalent and distinct from classic MD. This study has 
demonstrated a new type of mitochondrial dysfunction that may 
be the result of redox abnormalities and could affect a significant 
number of children with ASD. In fact, the LCL subgroup with 
mitochondrial abnormalities represented 44% of the total AD LCLs 
examined. 

Limitations 

The number of subgroups that could be identified depended on 
the total number of LCLs examined, which was limited. Future 
studies will need to examine a larger number of LCLs to confirm 
these findings and determine whether there are multiple LCL 
subgroups. In addition, studies need to verify these findings in 
better matched control samples. Age-matched controls LCLs will 
be very useful as they become available and matching to typically 
developing sibling LCLs may be an extremely powerful method to 
relate the mitochondria metabolism phenotype uncovered in this 
study specifically to autism. 



CONCLUSIONS 

This study has identified a novel pattern of mitochondrial 
dysfunction in immune cells derived from AD children that 
appears to be present in a significant subgroup of LCLs. Thus, we 
demonstrate a new type of mitochondrial disorder that may affect 
a significant subgroup of AD children and provide insight into the 
interactions between systems that have been independently 
demonstrated to be abnormal in ASD.^^ This information provides 
insight into the pathophysiology associated with ASD and a 
pathway for designing medical treatments for ASD. 
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